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We present a model for the spontaneous onset of fast magnetic reconnection 
in a weakly collisional plasma, such as the solar corona. When a current layer of 
macroscopic width undergoes collisional (Sweet-Parker) reconnection, a narrow 
dissipation region forms around the X-line. This dissipation region naturally 
becomes narrower during the reconnection process as stronger magnetic fields 
are convected toward the X-line. When the dissipation region becomes thinner 
than the ion skin depth, resistive magnetohydrodynamics breaks down as the 

c/> ' Hall effect becomes important and the Sweet-Parker solution ceases to exist. A 

transition to collisionless (Hall) reconnection ensues, increasing the reconnection 

>> rate by many orders of magnitude in a very short time. Predictions of the model 

are consistent with constraints set by observations of solar flares. 

Subject headings: Sun: flares -- Sun: corona -- magnetic fields -- plasmas - 

Sun: coronal mass ejections 
O 
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§ ■ 1. INTRODUCTION 

O' 

Magnetic reconnection is the driver of explosions in the solar corona. The first self- 
consistent description of magnetic reconnection, the Sweet-Parker model (Sweet 1958; Parker 
1957), was immediately recognized as yielding energy release rates too slow to explain obser- 
vations. Petschek reconnection (Petschek 1964) and related models (Priest & Forbes 1986) 
can be much faster, but requires anomalous resistivity (Sato & Hayashi 1979; Biskamp 1986), 
a process which is not well understood. In the past 10-15 years, a new paradigm of colli- 
sionless (Hall) reconnection has emerged, in which non-magnetohydrodynamic terms make 
reconnection much faster (Birn et al. 2001), about six orders of magnitude faster for solar 
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flare parameters. Signatures of Hall reconnection have been observed in magnetospheric 
observations (Oieroset et al. 2001; Mozer et al. 2002) and laboratory experiments (Cothran 
et al. 2003; Ren et al. 2005). 

However, explaining observed reconnection rates is only part of the problem. One must 
also explain why reconnection is explosive: explaining how magnetic energy can accumulate 
without significant dissipation and identifying the mechanism triggering the onset of fast 
(Hall) reconnection to release the stored energy are long-standing problems. 

Recently, Cassak et al. (2005) showed that a catastrophic transition from Sweet-Parker 
to Hall reconnection occurs when the width S of the Sweet-Parker dissipation region falls 
below the ion skin depth di = c/u P i, where to P i = (Airne 2 jm^) 1 ! 2 is the ion plasma frequency 
and n is the plasma density. This is the length scale at which magnetohydrodynamics 
(MHD) breaks down and the Hall term in the generalized Ohm's law becomes important. 
For 5 < di, the Sweet-Parker solution ceases to exist, and an abrupt transition to Hall 
reconnection ensues. 

The catastrophic transition to Hall reconnection when 5 ~ di was demonstrated by 
(externally) decreasing the resistivity rj. However, the idea that a solar eruption is caused 
by a change in the resistivity by a large enough amount to cause a transition from colli- 
sional to collisionless reconnection requires ad hoc assumptions about the poorly understood 
energetics of the corona. In this letter, we suggest that the transition can occur as the 
result of the dissipation region thinning due to the convection of stronger magnetic fields 
into the dissipation region during slow Sweet-Parker reconnection. This is a generic process 
which is entirely self-driven: it relies on no external forcing or fine tuning of any parameters. 
To our knowledge, this is the first self-consistent model for the spontaneous onset of fast 
reconnection. 

To see why transitions to fast reconnection are generic to the reconnection process, 

consider a macroscopic current sheet with a small but non-negligible resistivity. Without 

small scale structure, the Hall term in Ohm's law is unimportant, so the system undergoes 

Sweet-Parker reconnection. Since the resistivity is small, the dissipation region is embedded 

within the macroscopic current sheet. The width 5 of the Sweet-Parker current layer is given 

by (Sweet 1958; Parker 1957) 
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where ca up = B up /(4imiin) ' is the Alfven speed evaluated using the magnetic field B up 
just upstream of the dissipation region and L is the macroscopic length of the Sweet-Parker 
current sheet. During reconnection, stronger magnetic field lines are convected into the 
dissipation region by the inflow, causing a decrease in the Sweet-Parker layer width. When 
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the dissipation region becomes thinner than di a transition ensues. 

In the following Section, we review the salient properties of Sweet-Parker and Hall 
reconnection. In § 3, we describe the numerical simulations and their results, and we discuss 
implications for the onset of solar flares in § 4. 



2. SWEET-PARKER AND HALL RECONNECTION 

In the Sweet-Parker model (Sweet 1958; Parker 1957), a steady state is reached when 
the convective inflow of magnetic field lines is balanced by diffusion of field lines towards the 
X-line, 

where v in is the inflow speed. From continuity, v in L ~ v out S, where v out is the outflow speed. 
Therefore, the inflow Alfven Mach number Ma = Vi n /cA U p (a measure of the reconnection 
rate) is given by the right hand side of equation (1) since v ou t ~ CAup- The extreme elongation 
of the dissipation region due to i] being very small for most plasmas of interest throttles 
Sweet-Parker reconnection. For solar flares, the Sweet-Parker prediction of Ma ~ 2 x 10~ 7 
is far smaller than inferred from observations, where we used typical parameters of n ~ 
10 9 cm~ 3 for the preflare density, Bo ~ 100 G for the preflare coronal magnetic field, L ~ 
10 4 km for a typical flux tube length, and a classical resistivity of rj ~ 2 x 10~ 16 sec based 
on a temperature of T ~ 10 6 K (Priest & Forbes 2002). 

The physics of Hall reconnection is fundamentally different from that of Sweet-Parker. 
The motion of ions decouples from that of the electrons and the magnetic field at a distance 
di from the X-line. The electrons remain frozen-in to the magnetic field down to the electron 
skin depth d e = c/u pe . Where the species are decoupled, the Hall term in Ohm's law 
introduces whistler and/or kinetic Alfven waves (depending on the plasma f3) into the system 
(Mandt et al. 1994; Rogers et al. 2001). Both waves are dispersive with io oc k 2 . The 
dispersive property of these waves causes the outflow jet from the X-line to open as discussed 
by Petschek (Rogers et al. 2001). In the absence of dispersive waves, reconnection is slow as 
in the Sweet-Parker model. 

Numerical simulations (Shay et al. 1999; Huba & Rudakov 2004; Shay et al. 2004) have 
shown that the inflow speed for steady-state Hall reconnection is 

V in ~ O.lCAup (3) 

(i.e., Ma ~ 0.1). This result has been found to be independent of electron mass (Shay & 
Drake 1998; Hesse et al. 1999), system size (Shay et al. 1999), and dissipation mechanism 
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(Birn et al. 2001). Thus, we expect a dramatic increase in the reconnection rate when a 
transition from Sweet-Parker to Hall reconnection occurs. 



3. NUMERICAL SIMULATIONS AND RESULTS 

We perform numerical simulations using the massively parallel compressible two-fluid 
code F3D (Shay et al. 2004) in a periodic two-dimensional domain. The initial equilib- 
rium magnetic field is a double current sheet configuration given by one period of a sine 
sheet, B x o(y) = Bq cos(27ry / L y ) , where L y is the size of the domain in the inflow di- 
rection, with pressure balance enforced by a non-uniform density profile, n(y) = n + 
(B 2 /8nT ) sm 2 (2iiy/L y ). Here, n is a constant corresponding to the density at the edge 
of the domain and To = B 2 /4irno is the temperature, assumed constant and uniform for 
simplicity. The initial density at the center of the current sheet is, therefore, 1.5no. We 
impose no initial guide field. Lengths are normalized to the ion skin depth cZ i0 based on the 
density Hq at the edge of the computational domain, not the center of the X-line, which we 
denote as dix- Magnetic field strengths, velocities, times, and resistivities are normalized to 
Bo, the Alfven speed cao based on Bq and no, the ion cyclotron time fi^ 1 = (eBo/mic)^ 1 , 
and 770 = incAodio/c 2 , respectively. 

The computational domain is of size L x x L y = A09.6d i0 x 204.8d i0 with a cell size of 
0.1 djo x O.ldiQ. There is no viscosity, but fourth order diffusion with coefficient 2 x 10 -5 is 
used in all of the equations to damp noise at the grid scale. An electron mass of m e = rrii/25 
is used. Although this value is unrealistic, the electron mass only controls dissipation at 
the electron scales which does not impact the rate of Hall reconnection. A small coherent 
perturbation B x = —(0.004B L y /27r)z x V[sm(2irx/L x ) sm 2 (2ny/L y )} is used to initiate 
reconnection. The resistivity is taken to be uniform. Simulations are performed with i] = 
0.0025?7o and 0.0090?7o, both of which exhibit transitions to fast reconnection. We present 
results from the 77 = 0.0025i] simulation, which was initialized from the 77 = 0.0090r7o 
simulation at t — 5.364 kfi^ 1 . Initializing the simulation in this way introduces transient 
behavior, but it dies away (by t ~ 11 kfi^ 1 ) before small scale dynamics become important. 

When the system is evolved in time, the Hall effect is initially very small because the 
width of the current layer L y /2 = 102.4djo is large compared to dix, so the system evolves 
essentially as it would in pure resistive MHD. A Sweet-Parker current layer develops, as we 
will demonstrate later. The ion and electron inflow velocities, measured as the maximum 
value of the inflow into the X-line for each species, are plotted as a function of time late in the 
simulation in Figure la. Up until t ~ 18 kfi^ 1 , the electrons and ions are coupled as expected 
in MHD. The inflow speed is very small, but is slowly rising due to a gradual increase in 
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the upstream magnetic field strength B up as stronger magnetic fields are convected into the 
dissipation region. Figure lb shows the slow increase in B up , measured just upstream of the 
current layer in the simulation. 

When the ions decouple from the electrons, the inflow speeds begin to increase dramati- 
cally and the system begins a transition to Hall reconnection. This transition initiates when 
the width of the current layer S falls below dix, as is shown in Figure lc. The thick solid line 
is dix as a function of time. After decoupling, one must distinguish between the electron 
and ion current sheet widths, which we denote as S e and Si, respectively. The solid line is 
S e , determined by the half width at half maximum of the total current layer. The dashed 
line is Si, determined by the greater of S e and the half width at half maximum of the total 
inflow current. The latter becomes non-zero where the electrons and ions decouple, and is 
therefore a measure of the edge of the ion dissipation region. One can see Si decreasing from 
large scales (larger than dix) as the upstream magnetic field increases, and the transition 
begins when it is of the order of dix- 

Finally, to verify that the system is undergoing Sweet-Parker before the transition, and 
Hall reconnection after, we must check the validity of the inflow speed predictions from 
equations (2) and (3). The thick solid line of Figure Id shows Vi n as a function of time. 
The dashed line is the Sweet-Parker prediction from equation (2) (vi n ~ r]/S in code units), 
while the thin solid line is the Hall reconnection prediction with a constant coefficient of 
0.17, which is of the order of ~ 0.10 as expected from equation (3). Clearly, up until 
about t ~ 18 kfi" 1 , there is excellent agreement with the Sweet-Parker result. A grayscale 
plot of the current layer during the Sweet-Parker phase (at t = 11 A kfi" 1 ) is shown in 
Figure 2a, showing the characteristic elongated dissipation region (similar to those observed 
with pure MHD simulations by Jemella et al. (2004)). After a relatively brief transition 
time lasting until t ~ 19.5 kfi^ 1 , the inflow speed is well modeled by the Hall prediction. 
A grayscale plot of the current layer during the Hall phase (at t = 19.6 kfi" 1 ) is shown in 
Figure 2b, showing the open outflow configuration characteristic of Hall reconnection. We 
observe a large enhancement of the quadrupolar structure in the out of plane magnetic field, 
a signature of Hall reconnection (Mandt et al. 1994). Cuts across the current sheet at the 
X-line normalized to its maximum value are plotted as the dashed line and dot-dashed lines 
in Figure 2c, showing that S e falls to d e = 0.2rfj during Hall reconnection, as is expected 
when electron inertia provides the dissipation. For comparison, the solid line is a cut across 
the initial equilibrium current sheet. 
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4. DISCUSSION 

The spontaneous onset model presented here provides a possible explanation of why 
reconnection sites in weakly collisional plasmas are apparently quiet for a long time as 
magnetic energy accumulates before a sudden onset of fast magnetic reconnection releases 
it. A rigorous comparison of this model with flare observations is challenging because the 
dissipation regions associated with the transition to fast reconnection are much narrower 
than can be resolved with satellite or ground-based observations. We can, however, compare 
some basic predictions with observations. 

First, are macroscopic current sheets in the corona wide compared to the ion skin depth? 
Using values of the plasma parameters in a solar flare from § 2 gives an ion skin depth of 
only di ~ 7 x 10 2 cm, far narrower than expected macroscopic current sheets in the corona. 
At present, current sheets in the corona are inaccessible to observations, though reasonable 
scales for the current sheet width W s may be 100-1,000 km. 

Second, taking rj as a given, what is the critical upstream magnetic field strength B* 
which would make the Sweet-Parker current layer width equal to dp. Setting S = di in 
equation (1), we find 

B* ~ V4nm t n (^l\ ~5G (4) 

using the values from § 2. This is accessible during reconnection in the corona. 

Third, what is the time scale for the quiet time r q , during which Sweet-Parker recon- 
nection could be active but magnetic energy could accumulate? Since the field is frozen-in 
outside of the dissipation region, it is the time it takes for a field of strength B* to be 

convected in by the inflow, 

d^ 






(5) 



where £ is the distance upstream from the X-line. This can be approximated using V{ n ~ 
{r\(? /AticaL) 1 ^ 2 and by assuming a linear profile in the magnetic field B = Bq^/W s in ca- 
Integrating from £ ~ B*W s /Bq to approximately zero gives 



r,~2.r. l /-^-^~(-^-)4xI0««, (6) 

''r]c z CAoB \ 100 km/ 

where cao is the Alfven speed based on Bq. The numerical factor is about 11 hours, which 
is a reasonable time scale for the accumulation of magnetic energy due to footpoint motion 
in the photosphere (Dahlburg et al. 2005). 

The time it takes for the transition from Sweet-Parker to Hall reconnection, corre- 
sponding to the time from onset until maximum flare signal, can be bounded above by the 
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convective time across the Sweet-Parker current sheet S/vi n , which at the transition time is 
the same as the resistive time across the layer (tjc^/AttS 2 )" 1 and the convective time along the 
layer L/v out . For our simulation, the resistive time is ~ 400f2~ 1 , which compares reasonably 
well with the observed time of the transition (see Figure Id). For solar flare parameters, 
the resistive time across the layer is approximately 28 sec, which is comparable to the onset 
times seen in flares (Priest & Forbes 2002). The predicted observable parameters are quite 
consistent with solar flare phenomena. 

The present simulations do not include the effect of an out of plane (guide) field, the more 
generic configuration for magnetic reconnection. It was conjectured (Cassak et al. 2005) that 
the transition to fast reconnection in the presence of a guide field is also catastrophic, but 
occurs when the width of the current layer reaches the ion Larmor radius p s = c s /VL ci) where 
c s is the ion sound speed, instead of the ion skin depth d{. This is because p s is the scale where 
dispersive (kinetic Alfven) waves become important in the presence of a guide field (Rogers 
et al. 2001). Interestingly, recent laboratory experiments at the Versatile Toroidal Facility 
(Egedal et al. 2000) have observed spontaneous reconnection, and preliminary diagnostics 
suggest that the width of the current layer at onset is very close to their value of the ion 
Larmor radius p s (Egedal, private communication). 

Finally, Longcope et al. (2005) recently observed an active region reconnect with a 
nearby flux loop as it emerged from the corona. A phase of slow reconnection was observed 
for ~ 24 hr, during which magnetic energy accumulated in the corona. This was followed by 
fast reconnection lasting ~ 3 hr. The onset was sudden with no visible trigger mechanism 
observed. The energy released during fast reconnection was shown to be comparable to 
the energy accumulated during slow reconnection. Based on parameters inferred from the 
observations (a loop voltage of 10 9 V, a current sheet depth of 2 x 10 5 km, a sheet length of 
L ~ 3 x 10 4 km, a sheet current of / ~ 1.34 x 10 11 A, and a density of n ~ 10 9 cm -3 ), the 
fast reconnection rate was M& ~ 0.05, based on a reconnection electric field of E ~ 5 V/m 
and a reconnecting magnetic field of Bq ~ 4 G, consistent with Hall reconnection. These 
observations provide solid evidence for the accumulation of magnetic energy during a slow 
reconnection phase followed by a spontaneous onset of fast reconnection, as proposed here. 

This work has been supported by NSF Grant No. PHY-0316197 and DOE Grant 
Nos. ER54197 and ER54784. Computations were carried out at the National Energy Re- 
search Scientific Computing Center. 
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Fig. 1. — Time dependence (in thousands of ion cyclotron times) of the (a) ion (solid) and 
electron (dashed) inflow velocities, (b) upstream magnetic field strength B up , (c) electron 
(thin solid) and ion (dashed) current sheet widths 5 e and Si and ion skin depth (thick solid) 
at the X-line dix, and (d) ion inflow velocity from the simulation (thick solid), with Sweet- 
Parker theory (dashed, from eq. [2]) and Hall theory (thin solid, from eq. [3] with 0.17 
replacing 0.10). 
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Fig. 2. — (Color online) Grayscale plot of the current sheet during (a) Sweet-Parker recon- 
nection (at t = 11.4 kfi^ 1 ) and (b) Hall reconnection (at t = 19.6 kfi^ 1 ). (c) Cuts across 
the X-line for the same two sheets (dashed and dot-dashed, respectively) normalized to its 
maximum value. The initial current sheet profile is the solid line. Notice the color table for 
b) has been skewed for greater contrast and the amplitude of the current density is vastly 
different for the two sheets. 



